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Directional  Distributions  and  Mean  Square  Slopes  of  Surface  Waves 

Paul  A.  Hwang  and  David  W.  Wang 

Oceanography  Division,  Naval  Research  Laboratory ,  Stennis  Space  Center,  Mississippi 

Field  observations  show  that  the  crosswind  component  constitutes  a  significant  portion 
of  the  ocean  surface  mean  square  slope.  The  average  ratio  between  the  crosswind  and 
upwind  mean  square  slope  components  is  0.88  in  slick-covered  ocean  surfaces.  This  large 
crosswind  slope  component  cannot  be  explained  satisfactorily  based  on  our  present  mod¬ 
els  of  ocean  wave  directional  distributions.  Two-dimensional  spectral  analysis  of  3D 
ocean  surface  topography  reveals  that  bimodal  directional  distribution  is  a  common  fea¬ 
ture  for  wave  components  shorter  than  the  peak  wavelength.  The  calculated  result  of  the 
upwind  and  crosswind  mean  square  slope  components  using  a  bimodal  directional  distri¬ 
bution  function  is  in  very  good  agreement  with  field  measurements. 


1.  INTRODUCTION 

Ocean  surface  waves  are  the  roughness  element  of  the 
air-water  interface.  The  directional  distribution  of  the 
ocean  surface  roughness  is  an  important  parameter  in  air- 
sea  interaction  studies  and  ocean  remote  sensing  applica¬ 
tions.  Examples  include  magnitude  and  direction  of  energy 
and  momentum  fluxes  between  air  and  water,  scatterome- 
ter  wind  velocity  measurement,  synthetic  aperture  radar 
(SAR)  imaging  of  the  ocean  wave  field,  and  optical  detec¬ 
tion  of  surface  and  subsurface  features.  Close  correlation 
between  the  gas  transfer  velocity  and  surface  wave 
parameters,  especially  the  mean  square  slope,  has  been 
reported  extensively  in  the  literature  (e.g.,  papers  presented 
in  the  Air-Water  Gas  Transfer  Symposia,  Brutsaert  and 
Jirka ,  1984;  Wilhelms  and  Gulliver ,  1991;  Jahne  and 
Monahan ,  1995;  Donelan  et  al .,  2001). 

The  preferred  representation  of  the  ocean  surface  rough¬ 
ness  properties  is  the  directional  wavenumber  spectrum  of 
the  ocean  waves.  Such  data  become  available  only  re¬ 
cently.  The  results  on  the  roughness  wavenumber  spectrum 
come  in  small  pieces  because  of  difficulties  in  data  acquisi¬ 
tion  and  analysis.  Several  field  measurements  of 
wavenumber  spectra  of  short  waves  (from  several  millime¬ 
ters  to  several  decimeters)  obtained  by  scanning  laser  slope 
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sensing  have  been  reported  [e.g.,  Bock  and  Hara ,  1995; 
Hara  et  al,  1994,  1998;  Hwang  et  al ,  1996].  In  principle, 
scanning  slope  data  include  two  orthogonal  components  of 
ocean  surface  slopes  measured  in  the  space  domain  and 
detailed  information  on  the  directional  distributions  of  in¬ 
dividual  spectral  components  of  ocean  waves  are  available, 
but  only  limited  results  on  the  directional  properties  of 
short  waves  over  a  small  range  of  wind  speeds  have  been 
reported  [Hwang,  1995;  Klinke  and  Jahne ,  1995;  Hara  et 
al ,  1998]. 

Most  ocean  wave  directional  spectra  are  obtained  from 
temporal  measurements  acquired  by  directional  buoys 
[e.g.,  Longuet-Higgins  et  al ,  1963;  Mitsuyasu  et  al ,  1975; 
Hasselmann  et  al,  1973,  1980]  or  wave  gauge  arrays  [e.g., 
Donelan  et  al,  1985]  and  frequency  spectra  are  derived. 
The  transformation  from  frequency  to  wavenumber  domain 
is  carried  out  using  the  dispersion  relationship.  This  ap¬ 
proach  is  rather  successful  for  long  gravity  waves.  The 
application  to  shorter  waves  is  not  very  satisfactory  be¬ 
cause  of  the  large  Doppler  frequency  shift  caused  by  con¬ 
vection  of  short  waves  by  surface  currents.  The  large  fre¬ 
quency  shift  introduces  large  uncertainties  in  the  interpre¬ 
tation  of  the  length  scales  of  the  measured  (apparent)  wave 
frequency  [e.g.,  Phillips,  1985;  Hwang  et  al,  1996;  Done¬ 
lan  etal,  1999). 

For  surface  roughness  investigations,  the  surface  slope 
data  of  Cox  and  Munk  (1954)  remain  the  most  comprehen¬ 
sive  in  terms  of  the  range  of  wind  speeds  encountered  and 
the  extent  of  their  statistical  analysis.  Based  on  their  meas¬ 
urements,  the  ocean  surface  roughness  has  a  significant 
crosswind  component.  In  most  cases,  the  ratio  of  the 
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crosswind  and  upwind  mean  square  slope  components  is 
greater  than  0.7.  The  large  crosswind  to  upwind  ratio  is 
also  found  in  later  measurements  of  high  frequency  ocean 
wave  spectra  using  laser  slope  gauges  [e.g.,  Hughes  et  al , 
1977;  Tang  and  Shemdin ,  1983;  Hwang  and  S hem  din, 
1988].  This  ratio  is  much  larger  than  what's  expected  from 
applying  unimodal  directional  distribution  functions  cur¬ 
rently  established  in  the  ocean  wave  spectral  models. 

Banner  and  Young  [1994]  emphasize  that  bimodality  is  a 
robust  feature  of  the  directional  distribution  of  wind¬ 
generated  waves.  The  dynamic  process  that  produces  a 
bimodal  directional  distribution  is  the  nonlinear  wave- 
wave  interaction  mechanism.  Although  bimodal  directional 
distributions  have  been  obtained  from  directional  buoy 
data  using  analysis  techniques  such  as  maximum  entropy 
method  (MEM)  or  maximum  likelihood  method  (MLM) 
[Young,  1994;  Young  et  al,  \995;  Ewans,  1998;  Wang  and 
Hwang,  2000],  there  are  continuous  disagreements  on  the 
existence  of  directional  bimodality,  mainly  because  the 
quantitative  results  of  bimodal  directional  distribution  dif¬ 
fer  significantly  depending  on  the  processing  method  em¬ 
ployed.  Hwang  et  al  [2000a, b]  report  2D  wavenumber 
spectral  analysis  of  3D  surface  topography  measured  by  an 
airborne  topographic  mapper  (ATM,  an  airborne  scanning 
lidar  system).  The  directional  resolution  derived  from  3D 
topographic  data  is  excellent.  The  analysis  presented  in 
Hwang  et  al  [2000b]  illustrates  that  the  directional  resolu¬ 
tion  for  a  wave  field  generated  by  a  steady  10  m/s  wind 
field  is  better  than  10°  for  wavenumber  components  higher 
than  the  peak  wavenumber.  The  directional  spectra  of  the 
ATM  data  show  clear  directional  bimodality  [Hwang  et  al, 
2000b].  They  also  notice  that  cases  of  bimodal  directional 
distributions  are  commonly  found  in  the  directional 
wavenumber  spectra  derived  from  similar  3D  spatial 
measurements,  such  as  those  acquired  by  aerial  stereo  pho¬ 
tography  [Phillips,  1958;  Cote  et  al,  1960;  Holthuijsen, 
1983],  airborne  imaging  radar  {Jackson  et  al  1985]  and 
land-based  imaging  radar  [Wyatt,  1995]. 

In  this  paper,  we  present  computations  of  the  upwind  and 
crosswind  mean  square  slope  components  using  four  dif¬ 
ferent  directional  distribution  models.  These  calculations 
are  compared  with  the  measurement  of  mean  square  slopes 
by  Cox  and  Munk  [1954].  Due  to  the  fact  that  there  remain 
major  uncertainties  on  the  spectral  properties  in  the  short 
wave  regime,  we  limit  our  investigation  to  gravity  wave 
components  in  the  equilibrium  and  saturation  ranges  of  the 
surface  wave  spectrum.  Quantitative  results  of  the  ratio  of 
the  crosswind  and  upwind  slope  components  calculated 
using  four  directional  distribution  functions  are  presented 
in  Section  2.  Further  discussions  on  the  spectral  function 
and  directional  distributions  are  presented  in  Section  3.  A 
summary  of  the  study  is  given  in  Section  4. 


2.  DIRECTIONAL  DISTRIBUTIONS  AND  MEAN 
SQUARE  SLOPE  COMPONENTS 

2.  L  Field  Data 

The  mean  square  slope  results  reported  in  Cox  and  Munk 
[1954]  are  derived  from  analyzing  the  sun  glitter  patterns 
of  the  ocean  surface  obtained  from  an  aircraft.  The  area  of 
coverage  for  each  image  of  glitter  patterns  is  typically  on 
the  order  of  one-half  square  kilometer.  The  results,  there¬ 
fore,  yield  a  high  degree  of  statistical  confidence.  Based  on 
these  data,  the  total  mean  square  slopes  of  the  ocean  sur¬ 
face  increase  linearly  with  wind  speed,  and  the  following 
two  formulas  are  given: 

Sckan  =  5.12  X  1 0~3 1/  +  (3±4)xl0”3,  (1) 

and 

sxlick2  =  1 .56  X 1 0’3  U  +  (8  ±  4)  X 1 0"3 .  (2) 

In  (1)  and  (2),  the  wind  velocity,  U,  is  measured  at  12.5-m 
elevation.  Combining  their  results  with  other  later  datasets 
and  making  corrections  for  cases  that  are  affected  by  back¬ 
ground  swell,  a  slightly  different  formula  in  terms  of  U]0 , 
the  neutral  wind  speed  at  10-m  elevation,  is  given  in 
Hwang  [1997] 

sckJ  =  5 . 1 2  X 1  0"3  Ul0  + 1 .25  X 1  (T3 .  (3) 

The  data  of  special  interest  to  this  paper  are  those  collected 
from  slick  covered  surfaces.  Altogether,  Cox  and  Munk 
[1954]  report  9  slick  cases  (one  natural  slick  and  8  man¬ 
made  slicks)  with  wind  speeds  ranging  from  1.6  to  10.6 
m/s.  The  ratio  between  crosswind  and  upwind  components 
is  quite  large.  For  the  9  slick  cases,  the  mean  value  with 
one  standard  derivation  is  0.88±0.097,  and  none  of  the 
cases  has  a  ratio  less  than  0.75.  This  large  fraction  of  the 
crosswind  mean  square  slope  component  has  never  been 
explained  satisfactorily.  As  will  be  shown  below,  calcula¬ 
tions  using  established  unimodal  directional  distribution 
functions  under-predict  the  crosswind  component  consid¬ 
erably. 

2.2.  Directional  Distribution 

For  a  given  directional  spectrum,  %( k ,  &)=%(k)D(k,  O),  the 
upwind  and  crosswind  slope  spectra  can  be  expressed  as 

Xlu(k,e)  =  k2uX(k)D{k,e),  (4) 

and 

X\  c(k,d)  =  k2cx(k)D(k,d),  (5) 
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where  k  is  wavenumber,  subscripts  u  and  c  denote  the  up¬ 
wind  and  crosswind  components,  X\  W  is  the  slope  spec¬ 
trum,  relating  to  the  displacement  spectrum  %(k)  by 
Xxik)^]? %(k\  and  D(k,6)  is  the  directional  distribution 
function  (assuming  wind  direction  is  at  6  ~  0,  thus 
k=kco%6,  and  kc=ks\x\6).  In  the  following,  the  displacement 
spectral  function  is  assumed  to  be 


X(k)  = 


[ bu,g-°-5k -2'5, 

i  Bk~\ 


k  <  kf 

k  >  kj 


(6) 


where  £=4.6xl0'3,  6=5.2xl0"2,  and  kx  is  the  matching 
wavenumber  separating  the  equilibrium  and  the  saturation 
spectra;  the  magnitude  of  k}  is  in  the  neighborhood  of  6.5kp 
[e.g.,  Phillips ,  1977,  1985;  Hwang  et  al ,  2000a;  Hwang 
and  Wang  2000].  The  ratio  of  crosswind  and  upwind  mean 
square  slope  components  integrated  to  an  upper  limit 
wavenumber,  kh  can  be  written  as 


£  sin2  QD(k,0)d6 

dk 

jVrm 

f  cos2  6D(k,  9)d9 

J -n 

dk 

(7) 


For  individual  wavenumber  components,  the  ratio  is 


f  sin2  0D(k,6)d0 

m  =  ^~2 - •  (8) 

f  cos2  0D(k,0)d9 

y~iz 


Clearly,  the  specific  functional  form  of  the  directional  dis¬ 
tribution  plays  an  important  role  on  the  magnitudes  of  r(kf) 
and  d(k).  Four  different  directional  distribution  models 
reported  in  the  literature  are  investigated  here.  The  first 
three  models  [Mitsuyasu  et  al ,  1975;  Hasselmann  et  al , 
1980;  Donelan  et  al ,  1985,  with  modification  by  Banner , 
1990]  are  unimodal.  The  key  features  of  these  models  in¬ 
clude:  (i)  The  dominant  propagation  directions  of  all  wave 
components  are  in  the  wind  direction,  (ii)  The  directional 
beamwidth  is  narrowest  near  the  spectral  peak  and  be¬ 
comes  broader  as  the  wavelength  of  an  individual  spectral 
component  increases  or  decreases  from  the  peak  wave¬ 
length.  Of  these  three  unimodal  directional  distribution 
models,  the  Donelan  et  al  distribution  is  narrower,  and  the 
Hasselmann  et  al  distribution  is  broader,  than  the 
Mitsuyasu  et  al  distribution.  The  directional  distributions 
of  shorter  wave  components  (£>1.3^)  of  the  fourth  model 
[Hwang  et  al,  2000b,  Appendix  A]  are  bimodal.  As  stated 
earlier,  bimodality  is  a  robust  feature  of  nonlinear  wave- 
wave  interaction  [Banner  and  Young ,  1994].  The  beam- 
width  (an  integrated  property  of  the  directional  distribu¬ 


tion)  predicted  by  all  four  models  are  comparable  [Hwang 
et  al,  2000b]. 

Figures  la  and  b  show  the  computed  d(k)  and  r(k})  for  a 
wave  field  generated  by  10  m/s  wind.  The  ratio  between 
crosswind  and  upwind  slope  components,  d{k)  calculated 
using  the  bimodal  function  is  usually  larger  than  those  cal¬ 
culated  using  unimodal  directional  distributions,  especially 
in  the  lower  range  of  kJkp.  The  integrated  ratio,  r(kj),  based 
on  the  bimodal  function  is  much  closer  to  the  observed 
value.  Following  the  calculation  that  waves  shorter  then 
0.3  m  are  suppressed  by  slicks  [Cox  and  Munk,  1954;  Phil¬ 
lips,  1977],  the  corresponding  wavenumber  is  A/=27t/0.3 
rad/m.  For  this  wind  speed,  ks-2Hkp,  and  r(^>0.8.  In 
contrast,  r(ks)  is  approximately  0.7  based  on  the  models  of 
Mitsuyasu  et  al  [1975]  and  Hasselmann  et  al  [1980],  and 
is  less  than  0.6  based  on  the  model  of  Donelan  et  al 
[1985].  Hwang  and  Wang  [2000]  examine  the  effect  of 
cutoff  wavenumber  on  the  integrated  mean  square  slope, 
and  show  that  the  effect  is  minor.  Sensitivity  tests  on  other 
spectral  parameters  are  also  described  in  Hwang  and  Wang 
[2000]. 

The  ratio  r(ks )  as  a  function  of  wind  speed  is  plotted  in 
Fig.  lc.  For  comparison,  the  field  data  of  Cox  and  Munk 
[1954]  are  also  shown.  Of  the  four  directional  models,  the 
results  of  bimodal  directional  distribution  are  in  best 
agreement  with  the  field  data,  but  still  underestimate  the 
crosswind  to  upwind  ratio  by  approximately  14%,  the  other 
three  (unimodal)  models  underestimate  the  ratio  by  26  to 
39%. 


3.  DISCUSSIONS 

As  noted  in  Section  2,  the  crosswind  to  upwind  ratio  of 
mean  square  slopes  calculated  using  a  bimodal  directional 
distribution  is  still  14%  lower  than  the  field  data.  Several 
factors  may  contribute  to  the  observed  larger  crosswind 
surface  slope  component  in  the  ocean.  Firstly,  the 
computation  is  based  on  the  assumption  of  steady  wind 
forcing.  Fluctuations  in  the  wind  field,  especially  the  wind 
direction,  may  contribute  significantly  to  a  higher 
magnitude  of  r(ks).  Secondly,  the  derivation  of  the  bimodal 
directional  function  is  based  on  data  from  a  mature  wave 
field  where  the  dominant  wave  component  aligns  with  the 
wind  vector.  In  a  young  sea,  the  dominant  wave  may 
propagate  in  oblique  angles  with  respect  to  the  wind  vector 
in  order  for  the  wave  field  to  maintain  in  resonant 
condition  with  the  forcing  wind  field.  According  to 
Phillips  [1957]  resonance  mechanism  of  wind-wave 
generation,  two  wave  systems  propagating  at  oblique 
angles  symmetric  to  the  wind  are  generated.  Spatial 
measurements  of  3D  ocean  surface  topography  using 
airborne  scanning  radar  [Walsh  et  al,  1985,  1989]  and 
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Figure  1.  The  ratios  (a)  d(k)  (Eq.  8),  and  (b)  r(kj)  (Eq.  7),  calculated  for  a  wave  field  generated  by  10  m/s  wind,  (c)  The 
ratio  r(ks)  as  a  function  of  wind  speed,  where  the  upper  limit  of  integration  wavenumber,  ks=2n/03  rad/m,  representing 
the  wavenumber  above  which  the  waves  are  damped  by  the  slicks.  The  results  of  slick  cases  of  Cox  and  Munk  [1954] 
are  also  shown  as  triangles  in  the  plot.  Four  different  directional  distribution  models  arc  used  in  the  calculation:  dashed- 
and-dotted  curves:  Mitsuyasu  et  al.  [1975],  dotted  curves:  Hasselmann  et  al.  [1980],  dashed  curves:  Donelan  et  al. 
[1985],  and  solid  curves:  Hwang  et  al.  [2000b]. 


airborne  scanning  lidar  [Hwang  and  Wang ,  2000]  have 
obtained  examples  of  waves  propagating  in  oblique  angles 
with  respect  to  the  wind  vector  under  steady  wind  forcing. 
In  both  cases  presented  by  Walsh  et  al.  [1985,  1989],  most 
of  the  reported  directional  spectra  have  one  wind-generated 
wave  system.  At  shorter  fetches  (less  than  approximately 
100  km  for  10  m/s  wind  forcing),  the  direction  of  the  wind- 
wave  system  deviates  significantly  from  the  wind  vector. 
In  the  case  reported  by  Hwang  and  Wang  [2000],  two 
wave  systems  straddling  the  wind  vector  exist.  The  swell 
condition  in  the  dataset  of  Hwang  and  Wang  [2000]  is 
much  milder  than  those  encountered  in  Walsh  et  al  [1985, 
1989],  and  may  contribute  to  the  observed  differences.  In 
any  case,  one  expects  a  significant  increase  in  the  ratio 
between  crosswind  and  upwind  mean  square  slope 
components  when  waves  are  propagating  at  an  oblique 
angle  with  respect  to  wind. 


4.  SUMMARY 

The  mean  square  slope  dataset  of  Cox  and  Munk  [1954] 
is  regarded  a  masterpiece  in  the  study  of  ocean  surface 
roughness.  It  has  been  used  as  a  major  calibration  reference 
for  many  areas  of  research,  ranging  from  air-sea  interaction 
and  wave  dynamics,  to  acoustic  and  electromagnetic  re¬ 
mote  sensing  applications.  One  aspect  of  the  dataset  re¬ 
mains  quite  puzzling  over  the  last  half  century  or  so  is  the 
large  value  of  the  crosswind  slope  component.  For  slick 
covered  cases,  the  ratio  between  crosswind  and  upwind 
mean  square  slope  components  ranges  from  0.75  to  1.03. 
The  mean  value  with  one  standard  deviation  is  0.88±0.097. 
The  large  crosswind  slope  component  cannot  be  explained 
by  the  unimodal  directional  functions  established  in  mod¬ 
em  wave  spectral  models  (Fig.  lc).  The  mean  ratio  from 
three  unimodal  directional  distribution  functions  [ Donelan 
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et  al ,  1985,  with  modification  by  Banner,  1990;  Hassel- 
mann  et  al. ,  1980;  Mitsuyasu  et  al,  1975]  are  0.54+0.11, 
0.65±0.12  and  0.65±0.13,  respectively. 

Directional  spectral  analysis  of  3D  ocean  surface  topog¬ 
raphy  collected  by  an  airborne  scanning  lidar  system  shows 
that  in  a  mature  sea,  unimodal  directional  distribution  ex¬ 
ists  only  in  a  narrow  wavenumber  range  near  the  spectral 
peak.  For  wave  components  shorter  than  the  dominant 
wavelength,  bimodality  is  a  robust  feature  of  the  direc¬ 
tional  distribution.  Nonlinear  wave-wave  interaction  is  the 
mechanism  that  generates  the  bimodal  feature  [Banner  and 
Young,  1994;  Hwang  and  Wang,  2000].  Using  the  bimodal 
directional  distribution  function,  the  calculated  average 
ratio  between  crosswind  and  upwind  slope  components  is 
0.76±0.083,  which  is  in  much  better  agreement  with  field 
measurements.  This  mean  value  is  still  about  14%  lower 
than  the  field  data.  Possible  factors  contributing  to  the  ob¬ 
served  larger  crosswind  slope  components  in  the  field  data 
include  fluctuation  in  the  wind  field  and  less-than-mature 
stage  of  the  wave  field.  In  both  situations,  wave  compo¬ 
nents  traveling  at  oblique  angles  from  the  wind  vector  con¬ 
tribute  to  the  observed  larger  magnitude  of  the  average 
crosswind  slope  component. 

APPENDIX.  BIMODAL  DIRECTIONAL 
DISTRIBUTION  FUNCTION 


Hwang  et  al.  (2000b)  acquire  3D  ocean  surface  topogra¬ 
phy  using  an  airborne  scanning  lidar  system.  Based  on  data 
obtained  under  a  quasi-steady  wind  field  the  directional 
distribution  function  at  each  wavenumber  of  the  measured 
2D  spectrum  is  expressed  in  Fourier  series, 


D(k,0)  =  - 

71 


1  +  Z  A  (k) cosin 6 

n~  1 


-nl2<9<n!2' 


(Al) 


Coefficients  for  the  third  order  polynomial  fitting  to  each 
of  the  first  9  Fourier  components,  An,  n-\,  2,  9,  are 


tabulated  for  reconstructing  the  bimodal  function.  The  da¬ 
tabase  that  establishes  the  polynomial  coefficients  is  lim¬ 
ited  to  k  <  10  kp  so  extrapolation  too  far  beyond  10  kp  may 
produce  large  excursions  in  the  directional  function.  The 
computation  presented  in  this  paper  is  based  on  extrapola¬ 
tion  of  the  directional  coefficients  to  12 kp,  and  the  coeffi¬ 
cients  at  12 kp  are  used  for  the  remaining  higher  wavenum¬ 
ber  components  in  the  computation  range.  The  coefficients 
of  the  third  order  polynomial  fitting  of  the  first  9  Fourier 
components  are  tabulated  below. 
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